Draft version September 5, 2012 

Preprint typeset using I^-T^X style cmulatcapj v. 5/2/11 



(N 

o 

(N 
C/3 



CO 
Oh' 



> 
oo 

(N 

q; 

ON' 

o: 



X 



THERMAL TO NONTHERMAL ENERGY PARTITION AT THE EARLY RISE PHASE OF SOLAR FLARES 

Alexander A. Altyntsev 1 , Gregory D. Fleishman 2,3 , Sergey V. Lesovoi 1 , Natalia S. Meshalkina 1 

Draft version September 5, 2012 

ABSTRACT 

In some flares the thermal component appears much earlier than the nonthermal component in X-ray 
range. Using sensitive microwave observations we revisit this finding made by Battaglia et al. (2009) 
based on thorough analysis of RHESSI data. We have found that nonthermal microwave emission 
produced by accelerated electrons with energy of at least several hundred keV, appears as early as 
the thermal soft X-ray emission indicative that the electron acceleration takes place at the very early 
flare phase. The non-detection of the hard X-rays at that early stage of the flares is, thus, an artifact 
of a limited RHESSI sensitivity. In all considered events, the microwave emission intensity increases 
at the early flare phase. We found that either thermal or nonthermal gyrosynchrotron emission can 
dominate the low-frequency (optically thick) part of the microwave spectrum below the spectral peak 
occurring at 3-10 GHz. In contrast, the high-frequency optically thin part of the spectrum is always 
formed by the nonthermal, accelerated electron component, whose power-law energy spectrum can 
extend up to a few MeV at this early flare stage. This means that even though the total number of 
accelerated electrons is small at this stage, their nonthermal spectrum is fully developed. This implies 
that an acceleration process of available seed particles is fully operational. While, creation of this 
seed population (the process commonly called 'injection' of the particles from the thermal pool into 
acceleration process) has a rather low efficiency at this stage, although, the plasma heating efficiency 
is high. This imbalance between the heating and acceleration (in favor of the heating) is difficult to 
reconcile within most of available flare energization models. Being reminiscent of the trade off between 
the Joule heating and runaway electron acceleration, it puts additional constraints on the electron 
injection into the acceleration process. As a byproduct of this study, we demonstrate that for those 
cases when the optically thick part of the radio spectrum is dominated by the thermal contribution, 
the microwave spectral data yields reliable estimates of the magnetic field and source area at the early 
flare phase. 

Subject headings: Sun: filaments — Sun: flares — Sun: radio radiation — Sun: UV radiation 



1. INTRODUCTION 

The most prominent flare manifestations are explosive 
plasma heating and particle acceleration. The balance 
between thermal and nonthermal energy content varies 
during the flare and from one flare to another. The 
particle acceleration is widely believed to dominate at 
the impulsive phase; these electrons, after being accel- 
erated, propagate along closed magnetic field lines from 
coronal acceleration sites towards denser low corona and 
chromosphere. A direct diagnostics of the nonthermal 
electrons comes from hard X-ray observations of the 
bremsstrahlung as the flare-accelerated electrons collide 
with ions. The bremsstrahlung intensity is proportional 
to the ambient density; so the bulk of the hard X-ray 
(HXR) emission comes typically from the chromospheric 
footpoints where the precipitating flare-accelerated elec- 
trons lose most of their energy. In contrast, coronal 
HXR emission from the very acceleration sites is diffi- 
cult to detect except cases of acceleration in ver y dense 
flaring loops (jBastian et al.ll2007t IXu et alJl2008t ). How- 
ever, it is difficult here to isolate the directly accelerated 
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component because of essential in-situ Coulomb losses, 
whose effect depends on both transport regime and par- 
ticle energy. Another class of events in which the coro- 
nal HXR emission is easier to detect is the partially oc- 
cult ed_jlaxj5s_wfij3re^ are not visible 
(seelMasuda et al.lll995HWang et al.lll995l;lKrucker et al.l 
[2013)7 and a review bv lKrucker et al.l ()2008t) ). Analysis 
of such cases, however, suffers from the fact that the full 
flare energy balance is unavailable since the footpoint in- 
formation is missing. Fortunately, the situation is greatly 
improved if microwave observations are used to comple- 
ment the X-ray data. The radio emissions produced di- 
rectly at the coronal acceleration site can be firmly sepa- 
rated from other competing radio components, in favor- 
able cas es, allowing us to stud y the acceleration region 
directly (jFleishman et al.ll201lh . 

In addition to the main flare phase showing copious 
electron acceleration; observations provide strong evi- 
dence for a pre-flare coronal activity unrela ted to evap- 
orati on driven by precipitating electrons ([Harra et al.l 
|2001| ). Sometimes there is a prominent coronal soft X- 
ray (SXR) emission o ccurring several minu tes before the 
impulsive phase start (jFarmk fc S aw 199&f). Further im- 
portant knowledge about the pre-flare activity has been 
obtained from RHESSI spe ctral observations at t he low- 
energy range, 3-20 keV. iKrucker et all pOOl called 
this flare period with a pronounced gradual increase 
of the SXR emission occurring minutes before the de- 
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tectable onset of the HXR emission the early flare phase . 
The flares occu rred on 2002 July 23 (iLin et al.l 120031: 
Asai et alj|200ll and 2003 November 03 (jVeronig et al.l 
200l " represent vivid examples of flares with the early 
phase. 

The study of flare energization in the form of plasma 
heating and particle acceleration at this early phase is of 
particular interest, because it can clarify how the flare 
energy release is initiated and what mechanisms are re- 
sponsible for partition of the released energy between 
the thermal and nontherma l comp onents. For exam- 
ple, recently iBattaglia et al.l (|2009j ) have analyzed four 
flares and have demonstrated that morphology of these 
events showed one hot coronal source only, while the foot- 
points appeared much later, at the onset of the nonther- 
mal HXR emission. Remarkably, the X-ray spectra of 
these coronal sou rces at the early phase were nicely fit as 
purely thermal. IBattaglia ct al.l (|2009f) concluded that 
(i) the early flare phase is purely thermal; (ii) the de- 
tected increase of the thermal plasma emission measure 
(EM) requires some chromospheric evaporation, and (iii) 
the heat conduction is needed to drive this chromo s pheric 
evaporation. On the other hand Siarkow ski et al.l (|2009T j 
have studied the early phase in two simple flares, and 
concluded that the process of electron acceleration can 
occur during the early phase of flares; well before the 
impulsive stage. 

Microwave observati ons of the 2002 July 23 flare 
(jAsai et al.l [2001 [200l favor th e non-thermal alterna- 
tive (see, also. lWhite et al.ll2011l) . IHolman et all (|2003f ) 
examined HXR features of this flare, and reported that 
the nonthermal energy, even before the impulsive phase, 
was quite large. So, there is yet no clear and consis- 
tent picture of how the flarc/prcfiarc energy is released, 
and in what proportions it is divided between the ther- 
mal and nonthermal components. In particular, if the 
plasma is primarily heated by Coulomb losses of acceler- 
ated electrons, or directly by the energ y release remains 
an open, fundamental prob lem (see, e.g. lPetrosian fc Liu! 
[2001 iLhlfc Fletcherll2009( ). 

In this study we address this fundamental problem us- 
ing microwave data (from various instruments), which 
is more sensitive to the nonthermal electron component 
than the HXR data. Along with other available mi- 
crowave data, we use an essentially new input for the 
early phase study: the radio imaging at 5.7 GHz, which 
is indeed of primary importance since this frequency falls 
into the spectral peak region of the radio emission from 
the coronal sources at the early phase. Note, that the 
first observation of the 5.7 GHz early coronal source was 
performed back in 1992, well before the RHESSI era dur- 
ing X 9 limb flare on 1992 November 02 (jAltvntsev et al.1 
119991 ). The brightness center at 5.7 GHz was 10 Mm 
higher than the 17 GHz footpoint source. 

Observationally, this study can be facilitated by the 
fact that the preflare coronal sources are often reason- 
ably compact and weak. So the role of source inhomo- 
geneity and contribution from competing emission mech- 
anisms, otherwise essential, can be relatively minor in 
this case. Our event l ist co nsists of two (of four) events 
from IBattagli a et al.l (|2009l ) list complemented by three 
more events suitable for the analysis. No appropriate 
microwave data is ava ilable for two other events from 
IBattaglia et al.l (|2009l ) list. For the modeling of the 



microwave spectra we take into acco unt both thermal 
and nonthermal electron component s (|Benka fc Holmanl 
[19921: iFieisiiman fc Kuznetsovll20Tol) . 

In all the cases, we firmly detected a nonthermal elec- 
tron population as early as the thermal emission ap- 
peared. The energy content of the nonthermal electron 
population remained small and inessential for the plasma 
heating, indicative that the thermal plasma is heated by 
the flare energy release directly. In three (of five) cases, 
the low-frequency optically thick part of the microwave 
spectrum was entirely formed by SXR-producing thermal 
plasma, which offers a precise diagnostics of the source 
effective area and magnetic field. The source area (and 
thus, presumably the source volume) was found to grow 
in all these three events indicative that no chromospheric 
evaporation was needed to account for the EM increase. 
The observed EM increase is well consistent with a grow- 
ing source but without any change of the thermal elec- 
tron density. Although the overall nonthermal electron 
energy (and number density) content was small in all 
the cases, the acceleration mechanism was fully opera- 
tional, and powerful enough to produce a well-developed 
nonthermal power-law tails of the accelerated electron 
spectra, up to a few MeV, typical also for the impulsive 
flare phase. This implies that the key distinction between 
the 'normal' flare and the early flare phase is in remark- 
ably different efficiency of forming a seed suprathermal 
electron population via their extraction from the ther- 
mal pool. As soon as this seed population is formed, 
it is efficiently transformed to a power-law nonthermal 
distribution by the main acceleration process. Stated 
another way, the processes of electron injection from the 
thermal pool and their further acceleration are driven 
by distinctly different physical processes. The former 
process seems to be responsible for the energy partition 
between the thermal and nonthermal components, in the 
way that the stronger the plasma heating the weaker the 
injection, and vice versa. 

2. INSTRUMENTATION AND OBSERVATIONS 

The events under study were observed by a number of 
radio in struments: The Nobeyama Radio Polarimeters 
(NoRP, iTorii et al.l ll~979l ; iNakaiima efaLl [19851) record 
total and circularly polarized intensities (Stokes param- 
eters I and V) at 1, 2, 3.75, 9.4, 17, 35, and 80 GHz 
(Stokes I only). We also used the Radio Solar Telescope 
Network (RSTN) data. One second temporal resolution 
data were taken at 8 selected frequencies (0.245, 0.41, 
0.61, 1.415, 2.695, 4.995, 8.8, 15.4 GHz) from Learmonth 
station, Australia. 

The Nobeyama Radiohcliograph (NoRH, 

INakaiima et al.1 119941 ) obtains images of the Sun at 
17 GHz (Stokes I and V) and 34 GHz (Stokes I). The 
NoRH angular resolution is around 10 arcseconds at 17 
GHz, and 6 arcsec onds at 34 GHz Siberian Solar Radio 
Teles cope (SSRT; iSmolkov et ail [l98l iGrechnev et~aH 
|2003j ) produces microwave images at 5.7 GHz. In 
some events we also used one-dimensional brightness 
distributions measured with NS and EW arms of the 
SSRT. To analyze the most recent event in our event list 
(2011 February 15) in greater detail, we take advantage 
of more complete microwave observations obtained from 
new instruments developed by the SSRT team: the 
spectropolarimeter recording the total power data at 
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13 November 2003 




04:55:00 04:55:50 04:56:40 04:57:30 04:58:20 04:59:10 05:00:00 

Fig. 1. — Top: RHESSI corrected count rate light curves in 6- 12 keV (dashed), and 25-50 keV (dotted), during the early phase of the 
two events; same background is used as in IBattaglia eTaLl (120091) . Below: Microwave profiles recorded by NoRP, NoRH, and SSRT. The 
dotted vertical lines indicate the start and end of the early phase, the dashed lines show the time of the X-ray emission appearance at 
footpoints. 



2.3, 3.0, 4.0, 5.0, 6.0, 7.0 and 8.0 GHz (SP 2-9), and the 
10-antenna prototype of the multifrcquency heliograph 
mSSRT) observing at 5.1, 5.5, 6 3, 6 9 and 7.6 GHz 
Lesovov etaUl2009l: ILesovoi et al.M2012f ). 

We also use magnetograms from th e Michclson Doppler 
Imager (MDI, iScherrer etHI Il995[ ). and EUV images 
from the Extreme-Ultravio let Imaging Telescope (EIT, 
IDelaboudiniere et al.l 119951) on SOHO; as well as X-ray 
data from the Ra maty High-Energy Solar Sp ectroscopic 
Imager (RHESSI. IHn fe RHESSI Team|[200l and Fermi 
Gamma-ray Space Telescope (jMeegan et al.ll2009n . 

Flares on 2003 November 13 (N01E80) and 
20 04 January 6 (N05 E89). These two flares from 
the IBattaglia et afl (|2009| ) event list occurred during the 
SSRT observing time. The RHESSI corrected count rate 
light curves in 6-12 keV and 25-50 keV are shown in 
Figure [1] The early phase is marked by the dotted verti- 
cal lines (04:55:40-04:57:12 UT and 06:13-06:16 UT) and 
corresponds to the time interval of gradual increase of the 
6-12 keV light curve before appearance of the 25-50 keV 
emiss ion enhancement (see Figure 1 from[Battagli a~et al.l 
|2009|) . The microwave light curves recorded with the 
NoRP, NoRH and the SSRT are also shown. Note that 
the microwave emission at frequencies below 10 GHz ap- 
peared well before the impulsive phase. 

Images of the microwave sources at 5.7 GHz are shown 
in Figure [2] before the flare and at the early flare phase. 
Strong background microwave (presumably, gyroreso- 
nance) sources observed before the flares were situated 
close to the limb above AR 10537 and AR 10591, respec- 
tively (the dashed contours in Figure^]) . The correspond- 
ing microwave sources at the end of the early phase (the 
solid contours obtained after the background subtrac- 
tion) were located above the limb at di stances exceed- 
ing the half widths of the SSRT beam. IBattagli a et al.l 
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Fig. 2. — Contours: microwave sources before flares (dashed) and 
at the end of the early phase (solid). Contour levels are 9 X 10 5 K 
and 1.6 X 10 6 K for 13 November, and 3.4 X 10 5 K and 6.2 X 10 5 K 
for 6 January. The half widths of the SSRT beams, whose sizes 
and directions are dependent on the local time of observation, are 
shown by crossing bars in the top left corner. The half widths of the 
SSRT beam: 17 X 50 arcseconds (13 Nov) and 18 X 50 arcseconds 
(6 Jan). Line directions correspond to NS/EW scanning directions 
at times of the observations. The asterisks show centroid positions 
of the coronal X-ray sources from IBattaglia et al.l pOOSh . 

(|2009t) have studied the X -ray sources at this early phase 
in detail. Centroids of the X-ray loop top sources shown 
by the asterisks in Figure [2] are located clearly above 
the limb, implying the corresponding coronal loop top 
sources. The apparent sizes of the coronal microwave 
sources are larger than the corresponding X-ray sources 
because of relatively low spatial resolution of the SSRT, 
shown by the NS/EW crosses. The X-ray source volumes 
determined from the SXR imaging were 4.8 x 10 26 cm 3 
and 2.8 x 10 27 cm 3 by the end of the early stage, when 
their apparent sizes were 10.8" and 19.4", accordingly. 
During the early stage the source size remained roughly 
constant during 13 November early phase event, while 
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13 November 2003 06 January 2004 
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Fig. 3. — Microwave spectra at the end of the early phase and 
results of fitting by gyrosynchrotron emission generated by thermal 
and nonthermal electron distribution. Dashed curves show purely 
thermal gyrosynchrotron emission. Observations from NoRP (as- 
terisks) and SSRT (triangles). 

the area of the 06 January source increased, roughly, by 
a factor of three. The emission measures and temper- 
atures determined using the RHESSI observations were 
gradually growing up to 1.7 x 10 46 cm~ 3 and 25 MK for 
13 November event, and 1.1 x 10 46 cm~ 3 and 35 MK for 
06 January event. 

The microwave spectra observed at the end of the early 
phase are shown in Figure [3] The turn-over frequencies 
in both spectra were below the SSRT operating frequency 
implying that SSRT observed in the optically thin part 
of the microwave spectrum; i.e., the optically thin radio 
emission was produced from a coronal, not the footpoint, 
part of the source. Thus, wc can confidently conclude 
that the total power spectra obtained at other frequen- 
cies without spatial resolution pertain to coronal sources 
corresponding to the SSRT sources observed above the 
limb. 

X4.8 Flare on 2002 July 23 (S13E72). The coro- 
nal X-ray source appeared at the initial stage of this 
flare. It is the first proposed example of the early phase 
sources (|Krucker et all l2008f ). The early phase (00:18 
- 00:26 UT) continued until the HXR emission of the 
coronal source remained dominating. The X- ray obser- 
vation s at this phase were also a nalyzed bv | Lin et al.1 
([200l . IHolman et ahl (I2003T) . and ICaspi fc Linl (120101) . 
The coronal source displayed a gradual increase in both 
SXR and HXR emission below 60 keV. The RHESSI 
signal at hig h er en ergies appeared after 00:26 UT only. 
ICaspi fc Linl (|2010l ) have found three coronal electron 
populations: nonthermal, super-hot thermal, and hot 
thermal. Temperatures of the thermal electrons were 
determined at 00:20 UT and then were roughly constant 
and equal t o ^30 MK and ~ 20 MK, respectively (see 
Figure 1 in ICaspi fc Linl 120101 ) . The plasma density of 
the hot comp onent was es timated to be 1.7 x 10 10 cm -3 
(Fig.4.10 from lCaspill20T0t) . The HXR observations show 
that these plasma parameters did not change noticeably 
from 00:20 UT to 00:24 UT. However, the source volume 
was increasing from 3.7 x 10 26 cm -3 up to 7 x 10 26 cm~ 3 
during this same time interval. 

The microwave light curves (Figure 2|) grew gradu- 
ally during the early phase, and sharply afterward. The 
SSRT observations have began at the impulsive phase 
only, after 00:33 UT; although radio emission at frequen- 
cies up to 17 GHz appeared during the early phase. At 
the end of this phase the NoRP signal at 35 GHz ap- 
peared. A more sensitive Nobcyama Radioheliograph 
recorded the flare emission at 17 and 34 GHz well be- 
fore that time. After the 34 GHz emission appearance at 
00:22 UT the main microwave source was located slightly 
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Fig. 4. — Profiles of the integral microwave fluxes, recorded with 
the NoRP and RSTN. The dotted curves show the same fluxes 
divided by ten. The early phase bounds (00:18 - 00:26 UT) are 
marked by the dotted lines. 
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Fig. 5. — Spectra recorded with NoRP (asterisks) and RSTN 
(triangles) at the moments marked in Fig|4]by vertical dashed lines. 
The curves present results of spectral fitting by the TNT (thick) 
and the THM (thin) codes. 

above the SXR flare loop (|White et al.ll2003t lAsai et al.l 
120061 120091 ). The footpoint emission also appeared an d 
became essential after 00:22:30 UT (|Asai et al.l 120061 ). 
The spectrum evolution at the early phase is shown in 
Figure [5] Good correspondence between the NoRP and 
the RSTN data gives an idea on the measurement accu- 
racy. 

X7.1 Flare on 2005 January 20 (S06W90). Many 
studies were devoted to this extreme event, but they are 
mainly concerned the origin of accelerated protons, rel- 
ativistic electrons, an d the cosmic ray aspect (see, e.g., 
iGrechnev et al.l 120081 ). A gradual increase of SXR flux 
started at 06:00 UT, according to GOES-10 and GOES- 
12. At least half an hour before the flare maximum at 
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Fig. 6. — a) X-ray and microwave profiles in the early phase be- 
fore the flash phase started at 06:35:40 UT (marked by dash-dotted 
line). Dotted line illustrates the end of the Nobeyama observations. 
Top: GOES temperature and emission measure. Background mi- 
crowave and X-ray fluxes are subtracted. Bottom: NoRP and 
RHESSI light curves (in arbitrary units) in linear and logarithm 
scales, b) NoRP microwave profiles in the early phase. Bottom 
panel: natural logarithms of flux ratios in microwaves and X-rays. 
Vertical lines correspond to 06:20, 06:30:20, and 06:35:40 UT. 

06:46 UT, RHESSI detected the blob at loop topQ. 

The early phase lasted until 06:35:40 UT, where the 
temperature determined from GOES 10 data increases 
gradually from 7 to 13 MK without a noticeable increase 
of the emission measure (Figure . The HXR emission 
flux exceeds the background level at energies below 40 
keV, where it can be fitted by a power-law with index 
7.3 at 06:30 UT. 

First signatures of the nonthermal electrons appeared 
in microwaves after 06:20 UT. The radio fluxes at 2 - 17 
GHz were synchronously increasing during fifteen min- 
utes and reached a level of 100 SFU at 3 - 5 GHz. The 
microwave profiles are similar to X-ray profiles at ener- 
gies below 25 keV (see Figure |Bb, bottom). The degree 
of polarization of the total power emission was 10-20% 
LCP and almost did not change during the early phase 
at frequencies above 3.75 GHz. At 2 GHz the degree 
of polarization reached 30% RCP at some time frames. 
The similarity of the microwave and X-ray light curves 
favors a common origin of radio- and X-ray-emitting 
electrons. Sequence of the RHESSI images in the 3 - 
12 keV channels shows that the coronal source was lo- 
cated within the flare loop visible later in the impulsive 
stage (Figure [7]), and it was purely coronal until 06:38 
UT. The source volume was inc reasing during t h e first 
ten minutes of the early phase. IGrechnev et al.l (|2008f ) 
have estimated the loop parameters determined from 
EIT/171A and GOES/SXI images at 06:42:20 UT. The 
loop height was 4.6 x 10 9 cm, and width 9 x 10 s cm (vol- 
ume 7.8 x 10 27 cm 3 ). Plasma density was estimated to 
be (1 - 2) x 10 10 cm- 3 . 

Microwave sources are shown in Figure [7] at the end 
of the NoRH observations (06:30 UT). The 17 GHz and 
34 GHz emission sources are seen to originate at loop 
footpoints, located in the regions with opposite polar- 

4 |http://svs.gsfc.nasa.gov/vis/a000000/a003100/a003162/| 
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Fig. 7. — Flare structure in the early phase. MDI magnctogram 
(background, 05:46:30 UT); 3-12 keV RHESSI images (violet 
dash-dot-dotted contours at 0.3, 0.65, 0.8 levels, 06:41:00 - 06:41:30 
UT); 12 - 25 keV(turquoise dashed at 0.5, 0.65, 0.8 levels, 06:29:00 
- 06:32:00 UT); 5.7 GHz SSRT image (white solid contours at 0.5, 
0.9 levels, 06:33:26 UT); 17 GHz NoRH image in intensity (green 
continuous contours at 0.5, 0.7, 0.9 levels, 06:30 UT) and in the left 
circular polarization (green dotted at 0.9, 0.5 levels, 06:30 UT); 34 
GHz NoRH image in intensity (pink dash-dotted contours at 0.35, 
0.7, 0.9 level, 06:30 UT). The black dotted lines show the directions 
of the SSRT EW-array fan. The half widths of the SSRT beams 
are shown by crossing bars in the top left corner: 18 X 45". 
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Fig. 8. — One-dimensional distributions of the flare region at two 
moments during the early phase (solid and dotted curves). Vertical 
lines correspond to the EW fan directions in Fig[7| Lower solid and 
dashed curves show the subtraction of the two one-dimensional 
distributions in intensity and polarization, and correspond to the 
flare coronal source. Dashed curve shows distribution of circular 
polarization. 
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Fig. 9. — Microwave spectrum at 06:30 UT, recorded by 
NoRP (crosses), RSTN(asterisks), SSRT (squares), NoRH in the 
North footpont source (triangles), and the South source (dia- 
monds). The THT model spectra were calculated for loop top 
source (pink dashdot-dotted solid), the North footpoint source 
(green solid), and the South source (blue dash-dotted). The sum of 
the model spectra is shown by solid black thick curve. The dashed 
pink curve shows the contribution of gyrothermal emission from 
the coronal source. 

ities of the photosphere magnetic field. A line-of-sight 
magnetic field of f200 G (N-polarity) was measured in 
the South footpoint, while —1300 G (S polarity) in the 
North footpoint. At the early phase the flare emission 
at 5.7 GHz was relatively weak compared with the back- 
ground gyroresonant emission, which does not favor its 
direct imaging. However, the coronal source can be re- 
vealed by subtracting the background one-dimensional 
intensity distribution (Figure 18]). One-dimensional dis- 
tributions were recorded in the direction perpendicular 
to the fans of the linear arrays every 0.3 sec. The dotted 
lines in Figure [7] show the directions of the EW/SSRT 
array fan. 

Figure [8] presents one-dimensional distributions of flare 
region before and during the brightening at 06:30 UT. 
The subtraction of the former from the latter reveals the 
early flare phase source. The flare source was displaced 
to the West from the brightness center of the background 
source, so its location is likely to correspond to the loop- 
top region in Figure HI The FWHM size of the source is 
about 30", while the flux level of 28 sfu is about the total 
flare flux at this frequency. The degree of polarization 
reached 20% RCP; the brightness center of the polarized 
emission was displaced by 10" from the intensity bright- 
ness center to the North- West. Thus, the bulk of the 
5.7 GHz flare emission was produced from the coronal 
source. 

The microwave spectrum corresponding to the emis- 
sion enhancement at 06:30 UT is shown in Figure [9l 
The data at frequencies 2, 3.75, 9.4 GHz were recorded 
with the Nobeyama spectropolarimeter, at frequencies 
1.4, 2.65, 4.99, and 8.8 GHz with the Learmonth Solar 
Observatory (USAF RSTN Network) spectrometer. The 
SSRT flux at 5.7 GHz was determined from sequence of 
the one-dimensional EW distributions. It is reasonable 
to suppose that most of the emission below 6 GHz (see 
the radio spectrum in Figure [5]) originates from the coro- 
nal source. The microwave fluxes at 17 and 34 GHz, 
dominated by the footpoint sources (see Figure [7]) were 
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Fig. 10. — Top panel: RHESSI and Fermi X-ray signals at ener- 
gies 6-12 keV (pink dashed) and 25-50 keV(blue dotted). Back- 
ground microwave and X-ray fluxes are subtracted. Below: tempo- 
ral profiles recorded with spectropolarimeters NoRP, SP 2-9, and 
mSSRT. The vertical lines bound the early phase. 

obtained using the NoRH observations and are shown by 
triangles and diamonds in Figure [9] Assuming a power- 
law spectrum (F ~ J -7 ) at the range from 17 GHz to 
34 GHz, we can estimate ^North = 6.4 at the North foot- 
point, and Jsouth — 4-2 at the South one. 

X2.2 Flare on 2011 February 15 (N21W21). The 
flare on 2011 February 15 was the first X-class flare in 
the current solar activity cycle. Not surprisingly, it has 
already been analyzed in a number of papers. This flare 
occurred in AR 11158 in the central area of the solar disk. 
Given that the above study of the limb early flare phases 
have established that the main HXR (< 50 keV) and mi- 
crowave emissions (< 10 GHz) originate from compact 
coronal sources, it is reasonable to adopt the same mor- 
phology for the disk flares as well. 

For this flare analysis we take great advantage of new 
radio instruments developed by the SSRT team, and ob- 
serving now at the SSRT site. In fact, the 2011 February 
15 event is the first large flare recorded with these new 
SSRT instruments. The HXR and microwave profiles at 
the preflare phase are shown in Figure [TUJ Vertical lines 
bound an interval 01:45:10 - 01:47:55 UT. At the end of 
the early phase the estimates of plasma temperature and 
emission measure were 14 MK and 1.2 x 10 48 cm~ 3 from 
the GOES 15 X-ray data, and 34 MK and 3 x 10 46 c m - 3 
from the RHESSI data. The microwave emission was 
gradually rising during the early phase. Note that there 
is a good match between the light curves measured with 
different instruments. The emission at 17 GHz appeared 
during the early phase at 01:46:20 UT, while at 34 GHz 
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Fig. 11.— Flare on 2011 February 15 at 01:47:55. Microwave 
emission at 17 GHz: brightness temperature in intensity (back- 
ground), in circular polarization - white contours (right-handed - 
solid and left-hand - dotted). The HXR emission is shown in the 
6 — 12 keV by the pink dash-dotted contours, and in the 12-25 
keV by the green dashed contour. Axes are scaled in arcseconds. 
The half widths of the SSRT beams are shown by crossing bars in 
the bottom left corner: 19 X 25 arcseconds. 

even after this phase. 

The flare spatial relationships at the end of the early 
phase are shown in Figure 111! There was a loop seen 
in the 6-12 keV emission and in the microwaves. Its 
footpoints were clearly marked up by the oppositely po- 
larized sources at 17 GHz. The brightness center of the 
12-25 keV HXR source located at the loop top. The 
mSSRT observations have shown, through a forward fit 
of the measured visibilities, that the 4-9 GHz sources 
located near the loop top as well and their sizes were 
about 20" at the early phase. 

The microwave spectrum recorded at the end of the 
early phase with various instruments is shown in Figure 
[T2l There is a good correspondence between the differ- 
ent measurements, except some mismatch between the 
NoRH and RSTN fluxes at the frequencies above 15 GHz. 

3. MICROWAVE SPECTRUM MODELING 

A coronal source (e.g., a magnetic flux tube) can be de- 
tected via microwave observations in a few distinct cases: 
(i) via the thermal free-free emission (absorption) in case 
of high plasma density (large EM); (ii) via the thermal 
gyroresonant/gyrosynchrotron emission (GS) in case of 
high magnetic field and significant plasma heating; and 
(iii) via various nonthermal emission processes, of which 
we only consider below the incoherent gyrosynchrotron 
emission mechanism^ This means that having enhanced 
microwave emission does not necessarily imply that this 
emission is produced by a nonthermal electron popula- 
tion. In practice, the presence of the nonthermal compo- 
nent can be confirmed by microwave spectral modeling 
or forward fit, because the gyrosynchrotron spectra pro- 
duced by thermal or nonthermal electrons are distinctly 
different. Given that the early flare phase manifests itself 

5 Apparently, a combination of the mentioned mechanisms is 
also likely (see, e.e.. lBastian et al.ll2007l) . 
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Fig. 12. — Microwave spectrum at the end of the early phase 
(01:47:55 UT). 



as purely thermal in X-ray range ([Battaglia et al.ll2009D 
we checked if the microwave spectrum can be consistently 
described as being produced by the thermal component 
alone, and found that this purely thermal model fails in 
all the cases. This implies that a nonthermal electron 
population is always present. 

Therefore, it is reasonable to model the observed mi- 
crowave spectra by gyrosynchrotron emission generated 
by electrons, with a so-called thermal /nonthermal dis - 
tribution (TNT) proposed bv iBenka fc Holmanl (|1992l ). 
They demonstrated that the thermal plasma component 
can dominate the low-frequency optically thick part of 
the microwave spectrum. The TNT distribution 

!u THM (E), for E < E cr , 
A nth E~S, for E cr < E < E max . W 



ensures a smooth transition from the thermal to non- 
thermal distribution at a matching energy E cr . In the 
above expression, uthm(E) is the Maxwcllian distribu- 
tion function, and A nt h = UTHM{E cr )E s cr , to make the 
function continuous at the point E cr . In what follows we 
refer a nonthermal component to as all particles between 
E cr and E max . 

In practice we empl oyed fast gyrosynchrotro n code s 
recently developed by iFleishman fe Kuznetsovl (|2010f ) . 
which, by default, account for the corresponding free- 
free contribution along with the GS one. Given that 
available data is limited we are forced to adopt a few 
simplifications for our spectral modeling. In particular, 
we adopt the isotropic electron distribution in all the 
cases, and uniform source model in most cases (unless 
a different is explicitly stated). Thus, there are overall 
nine free parameters for the spectral modeling: T and 
nth, the temperature and the density of the background 
thermal plasma; matching energy E cr , power-law index 
5, the maximum nonthermal energy E max , the magnetic 
field strength B, the source depth D, and the area A; 
which characterize the emitting source, and the viewing 
angle 9 relative to the magnetic field dir ection. 

We start from two events from the Battagli a et al.l 
(|2009D list, whose morphology was studied in detail, 
so some guess parameters such as the source area A, 
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viewing angle, density and temperature of the back- 
groun d plasma can be tak en from the X-ray observa- 
tions (Battaglia et al. 2009). While the depths of the mi- 
crowave sources were estimated as \f~A~, uniform source 
model is expected to be valid for the compact looptop 
sources in these two flares. These input parameters ob- 
tained from the X-ray observations and used for the mi- 
crowave modeling are given by underlined italic charac- 
ters in Table 1. 

A classical GS spectrum displays a main spectral peak 
formed by the GS self-absorption at a frequency deter- 
mined by the condition r(/) ~ 1, where r(/) is the 
frequency-dependent optical depth of the GS source and 
can also contain a number of secondary (harmonic) peaks 
at small integer multiples of the gyrofrequency (see, e.g. , 
examples given in Figures by iBenka &; Holmanl 119921 : 
iFleishman fc Kuznetsovl 120101 and references therein) . 
The fast GS code is capable of reproducing the GS har- 
monic structure at relatively low frequencies; however, 
frequency gaps in the observations are too large to expect 
any reliable detection of the harmonic structure. In addi- 
tion, the harmonic structure is expected to be smoothed 
out by a relatively modest non-uniformity of the coro- 
nal source. T hus, we employ a continuous v ersion of the 
GS code ( see IFleishma n fc Kuznc tsov 20lT3, for the code 
nomenclature definition), which ensures high overall ac- 
curacy, while smoothies out the low-frequency harmonic 
structure. 

The results of the TNT fitting are shown by the thick 
curves in Figure while the dashed curves show the 
calculated spectra for the corresponding purely thermal 
electron distribution. Remarkably, the best fit spectra 
for these two events are distinctly different from each 
other. Indeed, for the 13 November 2003 event, the ther- 
mal contribution is negligible even if we adopt a relatively 
high magnetic held, 300 G, (a corresponding model spec- 
trum is shown in Figure GH left, dashed curve). This is 
in fact larger than needed to fit the high-frequency op- 
tically thin part of the microwave spectrum. The entire 
spectrum can be nicely fit by a nonthermal GS emission, 
if we adopt the source magnetic field to be above 70 G. 
The number density of nonthermal electrons n nt h was 
calculated using power law index S, for energies above 
E cr , and it is much less than nth- Nonthermal energy 
content was relatively small also (W n th/Wth "C 1). Note 
that the ratio of the plasma pressure to the magnetic 
pressure, commonly called the plasma (3, is about 0.2 for 
B = 70 G. 

In contrast, the thermal GS contribution (THM) en- 
tirely dominates the optically thick part of the microwave 
spectrum in the 06 January 2004 event; while the high- 
frequency optically thin part of the spectrum still has 
a clearly nonthermal origin. This thermal part of the 
GS spectrum is extremely conclusive in terms of the 
source parameter diagnostics. Indeed, in the optically 
thick regime, the brightness temperature of the thermal 
emission equals to the electron kinetic temperature T . 
Thus, the total power flux level is fully defined by the 
product A x T of the source area by the electron temper- 
ature, which immediately provides the effective source 
area A, given that the electron temperature is indepen- 
dently known from the SXR spectrum. Given that the 
plasma temperature does not change much during the 



studied early flare phase, being within 20-30 MK, the op- 
tically thick flux density increase can only be attributed 
to the corresponding increase of the effective source area 
and thus, the source volume. The volume increase can, 
in principle, be provided by either the source size or vol- 
umetric filling factor increase, from which we favor the 
former option given that the SXR imaging (see the pre- 
vious section) also reveals the SXR source size increase 
with time. In fact, the entire increase of the emission 
meas ure observed from the SXR data (Battagli a et al.l 
I2009T ) is accounted by this volume increase, and so no ad- 
ditional hot material (e.g., from a chromospheric evapo- 
ration) is needed. Then, for the given source size, density, 
and temperature, the spectrum peak frequency is solely 
defined by the source magnetic field. This offers simple 
and reliable diagnostics of the coronal source magnetic 
field. 

The 20 January 2005 event is somewhat different from 
the two events considered above, since its imaging data is 
inconsistent with the uniform coronal source model show- 
ing a strong footpoint contribution at high frequencies 
(17 and 34 GHz). This is in fact expected given relatively 
strong magnetic field values at the footpoints. Therefore, 
we explicitly consider three radio sources (one coronal 
and two footpoints), to roughly match the total power 
spectrum and also spatially resolved high-frequency data 
points. 

We modeled these two footpoint sources separately to 
match the corresponding 17 GHz and 34 GHz spatially 
resolved data. Recall, the microwave spectral indices 
are strongly different at the different footpoints, i.e., 
iNorth = 6.4 and t 'South = 4.2. The HXR footpoint spec- 
tral indices are also often different fr om each other, but 
the d ifference is typically within 0.6 (ISaint-Hilaire et ahl 
120081) . In contrast, the difference between the microwave 
indices is much larger, jNorth — Isouth = 2.2. A number 
of mechanisms has been proposed to account the HXR 
spectral index difference including anisotropic acceler- 
ation, different column densities in the opposite foot- 
points, differently nonuniform ionization, or effect of 
trapping due to m agnetic mirroring in an asymmetric 
magnetic trap (see ISaint-Hilaire et al.l [20081 for greater 
detail). Of which only the trapping can have an ef- 
fect on the microwave gyrosynchrotron spectrum. How- 
ever, formation of a noticeably asymmetric flaring loop 
is unlikely in our case, because the magnetic fields at 
the opposite footpoints arc comparable to each other 
in absolute value, being 1200 G and —1300 G, respec- 
tively. For the GS microwave spectrum from the foot- 
points, a much stronger effect on the spectrum slope is 
expec ted from the pitch-angle aniso tropy of the fast elec- 
trons (Fleishman & Melnikov 2003) having almost a pan- 
cake distribution at the mirror points. However, because 
we do not have a reliable ground to treat this anisotropy 
quatitatively, we assume the isotropic velocity distribu- 
tion of fast electrons, and match the high-frequency data 
points with Ssouth = 4.3 for the southern footpoint, while 
with &North — 5 for the northern footpoint. 

Given the above note on the possible effect of the pitch- 
angle anisotropy, the found mismatch between Ssouth and 
o~North does not necessarily imply a real difference in the 
fast electron distributions over energy in the different 
footpoints, but can simply indicate inadequacy of the 
adopted assumption of the distribution isotropy, which 
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is not too important for our study. However, this does 
imply that the fast electron energy spectrum can be 
harder than those derived from the footpoints, S < 4.3. 
Here, unlike other cases, we account the GS harmonic 
structure to give an idea of its possible role in the total 
spectrum formation. The thermal plasma contribution is 
negligible, since the fast electron number density is rela- 
tively large in this event. The input (underlined italic), 
and best-fit parameters are presented in Table [TJ 

Two other events display microwave spectra overall 
similar to that of the 06 January 2004 event, having ther- 
mal low-frequency parts and nonthermal high-frequency 
parts. In addition, we studied the spectrum evolution 
during the early stage of the 23 July 2002 event, and 
found that the thermal emission of an optically thick uni- 
form source well fits low frequency part of the observed 
spectra until 00:23:45 UT. Note that the assumption of 
the source homogeneity becomes incorrect at the end of 
the early phase because of essential emission c ontribution 
from footpoints seen in 17 and 34 GHz images (|Asai et al.l 
I2009T ). Perhaps it ca lls for three dimension al simulations 
of the flaring loops (jKuznetsov et al.ll2011| ). 

The spectral modeling of the 23 July 2002 yields the 
source magnetic field around 200 G. Interestingly, the 
flare magnetic field was independently estimated from 
polarizati on of the optic ally thin free-free microwave 
emission (|Asai et al.ll2009f) just before the flare. The de- 
gree of the circular polarization measured with NoRH at 
17 GHz was about 1.6%. The degree of polarization is 
linearly proportional to the line of sight magnetic field 
component. Thus, being combined with the determined 
magnetic field value of 200 G, it yields an estimate of 
the viewing angle at the source to be 8 = 75°, which 
is reasonable for a coron al so urce. From the RH ESSI 
observations. iCaspH (|2010T l and lCaspi fc Linl (|201Q[ ) have 
concluded that a hot thermal plasma with temperature 
about 30 MK, and density nth > 10 10 cm -3 , was present 
at the beginning of prc-impulsive phase in the corona. 
The parameters obtained from the microwave modeling 
for these events are given in Table 1 . The estimates from 
X-ray observations are available after 00:22 UT (they are 
given in brackets in Table 1). The density and tempera- 
ture estimates obtained from X-ray data agree well with 
those of the microwave spectrum modeling. The source 
area increase was detected in both X-ray and microwave 
domains, although the microwave source sizes were no- 
ticeably larger than the SXR source sizes. This is ex- 
pected because the SXR emission is optically thin and 
so the SXR source size represents, in fact, a FWHM of 
the emission measure distribution, while the microwave 
source area relates to optically thick emission and so rep- 
resents a full area of the source defined as a locus of all 
points corresponding to lines of sight with r > 1. 

Finally, the spectrum shape of the 15 February 2011 
event is similar to the 06 January 2004 and the 23 July 
2002 events. Although, it has a noticeably higher spec- 
trum peak frequency indicative of correspondingly larger 
magnetic field, ~ 500 G. For this event we used two input 
sets obtained, respectively, from the GOES and RHESSI 
(in brackets) data. Parameters resulted from the mi- 
crowave spectral modeling are summarized in Table [TJ 
including energy partition between the magnetic, ther- 
mal, and nonthermal energies. 



4. DISCUSSION 

So far, the simultaneous primary heating and acceler- 
ation are rarely studied together. It has been implicitly 
assumed that one of these processes is entirely dominant 
over the other. For example, it is commonly believed 
that the particle acceleration is a dominating process in 
a flash flare stage. Let us discuss now what microwave 
data says about the early flare phase in addition to the 
commonly used X-ray data. 

First of all, in all events under study there are non- 
thermal electrons that generate gyrosynchrotron emis- 
sion at frequencies above the spectral peaks. In the 
case of power-law distribution of emitting electrons the 
best fit indices are in the range from 2.5 up to 4, with 
the high-energy cutoff above 1 MeV. The densities and 
energy contents of nonthermal components (i.e., above 
E cr ~ 10 — 20 kev, see Table 1) were well below the ther- 
mal density and energy of the coronal sources. In two 
of the studied cases the emission of nonthermal electrons 
also dominates the observed spectra at frequencies below 
the spectrum peaks. 

Secondly, we found that thermal GS emission domi- 
nates the low-frequency microwave spectra in many cases. 
This offers reliable diagnostics of the source area and 
magnetic field. The radio estimate of the coronal source 
area is highly important because it is unbiased by the 
plasma density distribution, which is unlike the SXR- 
derived source area. Based on the radio data we confi- 
dently conclude that the source area grows at the course 
of the flare This fully accounts for the observed increase 
of the SXR-dcrived emission measure, while no density 
increase is needed. This means that no essential chro- 
mosphcric evaporation occurs in the analyzed cases, so 
no energy deposition to the chromosphere in the form 
of either precipitating electrons or heat conduction takes 
place. 

Thirdly, even though the thermal plasma contribution 
to the microwave spectrum is often essential, no purely 
thermal stage has been detected. Indeed, radio signa- 
tures of the nonthermal particles appear as soon as the 
plasma heating. Thus, the RHESSI (or Fermi) non- 
detection of the nonthermal emission at the early flare 
phase is accounted by its relatively low sensitivity, while 
the microwave observations turn out to be more sensi- 
tive to small numbers of the nonthermal electrons. Note, 
that in some cases (e.g., when the source area is some- 
what small, see the 13 November 2003 event as an exam- 
ple) , the entire microwave spectrum is dominated by the 
nonthermal GS contribution. 

The plasma beta /? is smaller than one, while the non- 
thermal energy density is much smaller than the thermal 
one in all the cases (i.e., W nt h *C Wth "C Wb)- Al- 
though the total energy content of the accelerated elec- 
trons is small, the available nonthermal electrons are high 
efficiently accelerated from slightly nonthermal to rela- 
tivistic energies. Their spectra are hard, 8 = 2.5 — 3 in 
most cases, and extended up to a few MeV. Stated an- 
other way, the shape of accelerated particle spectrum at 
the prcflarc phase is similar to that during flares, even 
though their levels (normalizations) are highly different. 
This can happen, for example, if the same, presumably 
stochastic, acceleration mechanism capable of accelerat- 
ing the charged particles from somehow created "seed 
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population" is involved at both preflare and flare phases. 
However, these seed populations must be formed differ- 
ently in prcflares or flares. 

The observed significant plasma heating suggests that 
the corresponding flare energy is already available. How- 
ever, it is divided highly unevenly between the plasma 
heating and nonthermal seed population creation. It 
works in a way similar to that in the presence of a 
DC electric field: at a preflare phase, a relatively large, 
but still essentially sub-Dreicer field, will heat the ambi- 
ent plasma via the quasi- Joule dissipation (an enhanced, 
anomalous resistivity is needed to yield a significant 
plasma heating), while the fraction of the runaway elec- 
trons capable of forming the mentioned seed population, 
will remain relatively minor. For a larger DC electric field 
the fraction of the runaway electrons will grow quickly, 
resulting in a powerful nonthermal component needed to 
produce the impulsive flare phase. 

Let us estimate what DC field E is required to form 
the seed populations at the preflare phase. Adopting 
typical parameters, n t h ~ 10 10 cm" 3 and T ~ 30 MK 
(see Table [T]), the electron Dreicer field is about E]j e m 
3 x 10~ 5 V/cm. The nonthermal to thermal electron 
number density ratio is about 10~ 4 (see Table [T]). To 
build this nonthermal component from the maxwellian 
tail, electrons with v > v cr (where exp(— v^ r /2v^ h ) ~ 
10~ 4 ) must runaway due to the DC electric field. Given 
that v 2 cr ~ (E De /E)v? h , we find E ~ 1.5 X 10~ 6 V/cm. 
Over a typical source size of ~ 10 9 cm, an electron can 
gain about 1 kcV of energy. This energy is far too small 
compared with the observed electron energies of 1 MeV 
or above. So this assumed DC field plays no role in form- 
ing the nonthermal power-law distribution responsible 
for nonthermal GS radiation. However, this ~ 1 keV of 
energy gain can be sufficient enough to form a slightly 
suprathermal seed population, from which the bulk (pre- 
sumably stochastic) acceleration produces the observed 
nonthermal power-law tails up to relativistic energies. 

Microwave observations, therefore, show that the en- 
ergy release mechanism in the preflare phase is accompa- 
nied by particle acceleration. Although this acceleration 
is much milder than that in the impulsive phase, and 
therefore, is hard to detect in X-rays. The nonthermal 
emission produced by accelerated electrons with energy 
of several hundred keV to a few MeV appears as early 
as the soft X-ray emission. Thus the non-detection of 
the hard X-rays at that early stage of the flares seems to 
originate from a limited sensitivity of available X-ray in- 
struments. The frequency of the spectrum peak is below 
10 GHz for the early flare phase of microwave emission in 
all cases, because of a relatively small number of accel- 
erated electrons at the radio sources. Microwave spectra 
show that magnetic field in the coronal sources are a few 
hundred gauss at the early phase. In some cases the num- 
ber of accelerated electrons is so low that the gyrosyn- 
chrotron emission from thermal electrons dominates the 
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low frequency part of the microwave spectrum. In these 
cases microwave spectrum shape provides the magnetic 
field estimate in the coronal sources. The microwave ob- 
servations of such events are promising for studying the 
transitions from the gradual energy release to the flash 
flare explosives. 

5. CONCLUSION 

The described findings give rise to a number of funda- 
mentally important conclusions about the flare heating 
and acceleration, (i) The flare energy release is capable of 
directly heating the thermal plasma, without noticeable 
in situ heating by fast electron beams or chromospheric 
evaporation driven by either electron beams or heat con- 
duction. Since no energy transfer process is detected, 
we propose the early flare phase sources to represent 
the energy release and acceleration sites, (ii) The effi- 
cient acceleration process of only a minor fraction of the 
plasma electrons implies that the acceleration mechanism 
involved is not capable of accelerating electrons directly 
from the thermal pool, but requires a pre-extracted (in- 
jected) seed electron population, (iii) This implies that 
the electron injection from the thermal pool and their 
further acceleration toward higher energies are driven by 
physically distinct processes. The first of them is ineffi- 
cient or somehow suppressed during the early flare phase, 
while the second is already fully operational. In par- 
ticular, the acceleration by cascading turbulence alone 
seems to be insufficient here. Since in the correspond- 
ing acceleration model both injection and acceleration 
are driven by the same turbulence intensity, so having 
inefficient injection while efficient acceleration looks at 
odds to this acceleration model, (iv) The observed sig- 
nificant plasma heating suggests that the corresponding 
flare energy is already available. However, it is divided 
highly unevenly between the plasma heating and non- 
thermal population creation. We propose, this is due to 
yet unspecified energy partition process operating in a 
way showing some resemblance to that controlling the 
balance between Joule heating, and runaway electrons in 
a DC electric field. We emphasize, that in addition to 
this energy partitioning process, some sort of stochastic 
acceleration capable of producing the observed power- 
law electron spectra is still needed during both early and 
impulsive flare phases. 
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TABLE 1 

Parameters of radio sources in five events. 
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